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Abstract 

We perform unbiased tests for the dumpiness of the Universe by confronting the Zel’dovich- 
Kantowski-Dyer-Roeder luminosity distance, which describes the effect of local inhomogeneities 
on the propagation of light with the observational one estimated from measurements of standard 
candles, i.e., type la supernovae (SNe la) and gamma-ray bursts (GRBs). Methodologically, we 
first determine the light-curve fitting parameters which account for distance estimation in SNe la 
observations and the luminosity/energy relations which are responsible for distance estimation of 
GRBs in the global fit to reconstruct the Hubble diagrams in the context of a clumpy Universe. 
Subsequently, these Hubble diagrams allow us to achieve unbiased constraints on the matter 
density parameter as well as the dumpiness parameter g which quantifies the fraction 
of homogeneously distributed matter within a given light cone. At a la confidence level, the 
constraints are Dm = 0.34 ± 0.02 and g = 1.00lg Q 2 from the joint analysis. The results suggest 
that the Universe full of Friedman-Lemaitre-Robertson-Walker fluid is favored by observations 
of standard candles with very high statistical significance. On the other hand, they may also 
indicate that the Zel’dovich-Kantowski-Dyer-Roeder approximation is a sufficiently accurate 
form to describe the effects of local homogeneity on the expanding Universe. 
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I. INTRODUCTION 


The standard physical model of cosmology is based on the solntion of general relativity 
describing a spatially homogeneous and isotropic spacetime, known as the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) solution. It is assumed that the geometry of our 
Universe is smooth on large scales. One of the major tasks in modern cosmology is to 
precisely determine the parameters which characterize the postulated model by htting 
the observational data. The cornerstone of observational evidence that supports the 
FLRW model is the existence of highly isotropic cosmic microwave background radiation 
(CMBR). It could be inferred that the spacetime should be exactly FLRW when the 


Q. 


background radiation appears to be exactly isotropic to a given family of observers 
Therefore,we can prove the Universe to be FLRW just from our own observations of the 
CMBR by taking the Copernican principle into consideration. Moreover, this result could 
be extended to the case of an almost isotropic background radiation, which hints at an 
almost FLRW spacetime j^. Although this simple solution of Einstein held equation 
provides an excellent description for the universe on large scales, it also makes clear 
that we need to understand the departures from a spatially homogeneous model when 
interpreting observational data. Indeed, departures from perfect homogeneity change the 
distance-redshift relation. However, in practice, cosmological observations are usually 
htted just using relationships derived from homogeneous models. 

The fact that matter is not continuously distributed can imprint most cosmological ob- 


S). 


servations probing quantities related to light propagation(as discussed in detail in Ref. 
in particular regarding the propagation of light with narrow beams, such as the redshift, 
the angular diameter distance, the luminosity distance, and the image distortion. The 
importance of quantifying the effects of inhomogeneities on light propagation was hrst 
pointed out by Zel’dovich jj] and Kantowski j^. They designed an “empty beam” ap¬ 
proximation by arguing that photons should mostly propagate in vacuum. Later, this 


was generalized by Dyer and Roeder as the “partially hlled beam” approach 6 


generally, the early work of Ref. j4| stimulated many studies on this issue 


0-2O|. 


7|. More 
In this 


framework, the proportion of clumped matter with respect to the homogeneous fluid is 
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characterized by the dumpiness or smoothness parameter. In addition, they arrived at an 
equation for the angular diameter distance which, via the Etherington relation, connects 
to the observable luminosity distance. We refer to it here as the Zel’dovich-Kantowski- 
Dyer-Roeder (ZKDR) luminosity distance. 

Since the 1960s, a rich literature has formed which concerns the ZKDR approach and its 
cosmological implications. Phenomenologies and investigations involving many different 
physical aspects were performed, such as analytical or approximate expressions 
critical redshift for the angular diameter distance 2^, gravitational lensing 


25 


21 


23|, 


accelerated expanding Universe models driven by particle creation 


26| . and 


27| . Recently, some 


quantitative analysis from such compact radio sources as standard rulers 28|, l29L and such 
type la supernovae (SNe la) or gamma-ray bursts (GRBs) as standard candles [30H35[ were 
also performed. To be specihc, in Ref. 30|, constraints on the dark energy and smoothness 
parameter from the so-called gold SN la sample released by the High-z Supernova team 


and the first year results of the Supernova Legacy Survey (SNLS), which is a planned 
five-year project [37|, were examined. The results suggested that SNe la data alone was 
incapable of constraining the smoothness parameter although the gold SN la provided a 
little more stringent constraint since this sample extended to appreciably higher redshifts. 
Later, Busti et al. 3^ performed an updated investigation where the statistical ana! 


was based on the 557 SNe la Union2 compilation data 


and 59 Hymnium GRBs 


ysi s 


39|, 


and almost the same conclusion was achieved. More recently, this issue was also studied 
by using Union2.1 SN la 40| plus nine long GRBs in 1.55 < ^ < 3.57 4l|] and the 
constrained value of the smoothness parameter indicated a clumped Universe 33|]. On 
the other hand, as concluded in their work, this result may be an indication that the 
ZKDR approximation is not a precise form of describing the effects of dumpiness in the 
expanding Universe. 

However, in these previous analysis, all distances of SNe la and GRBs applied to 
test the inhomogeneity of the Universe were derived from a global £t in the context 
of standard dark energy scenarios where the dumpiness has vanished, i.e., the flat A 
cold dark matter (ACDM) or tcCDM model. That is, the light-curve htting parameters 
accounting for the distance estimation in SNe la observations (e.g., a and /3 in the most 
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widely used SALT2 training method are left as free parameters (on the same weight 
as cosmological parameters) and are determined by htting the distances of SNe la, which 
is a linear combination of light-curve htting parameters and observed quantities, to the 
model-predicted ones in the context of the standard ACDM or rcCDM scenario. Therefore, 
HDs constructed in this way are somewhat model dependent. Moreover, cosmological 
implications on nonstandard dark energy scenarios or a Universe with homogeneity taken 


into consideration derived from these HDs are model biased 


43| . It has been shown 


that this kind of 
cosmology 


44 


oias cannot be neglected and may be signihcant in the era of precision 


45| . Certainly, this kind of bias also hides in the GRB cosmology where 


luminosity relations being responsible for distance es 


the mo del-dependent HDs of low-redshift SNe la 46 


imation of GRB are calibrated with 


47|. 


In this paper, we hrst reconstruct Hubble diagrams for the latest SNe la and for long 
GRB observations by calibrating the light-curve htting parameters and luminosity rela¬ 
tions, respectively, in the context of an inhomogeneous Universe with the cosmological 
constant. These Hubble diagrams can lead to unbiased tests for the matter density pa¬ 
rameter Vtm as well as the dumpiness parameter r]. For the joint light-curve analysis of 
the SDSS-H and the SNLS (JLA SN la) in the range of 0.01 < < 1.23 j^|, the con¬ 
straints are Dm = 0.291 o;q 5 and rj = slightly indicating a clumped Universe. For 

the long GRBs in the range of 1.48 < 2 ; < 8.20 49|, the best hts are = 0.42 ± 0.06 


and Tj = 1 . 00 ^q; 52 , strongly supporting a homogeneous Universe. For the combination of 
these two probes, the constraints are Dm = 0.34 ± 0.02 and r] = 1.001q;o 2, also favoring a 
universe full of FLRW huid with a very high conhdence level. We suggest that the matter 
density parameter Dm is mainly determined by the SNe la observations while the dumpi¬ 
ness parameter t] is primarily constrained from the observed GRB events. Moreover, it is 
also shown that larger scales are explored, the test more strongly implies a homogeneous 
Universe. These reasonable results may be an indication that the ZKDR approximation 
remains to be a precise description for the luminosity distance-redshift relation in a locally 
inhomogeneous Universe with the cosmological constant. 
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II. THE ZKDR LUMINOSITY DISTANCE 


For most cosmological models, angular or apparent size distance, which is proportional 
to the square root of the cross-sectional area A{z), is related to the luminosity distance 
by dA^z) = di^{z)/{l + z^. In the model only including dark matter and dark energy, 
the luminosity distance di,{z), which accounts for a partially depleted mass density in 
the observing beam but neglects lensing by external masses, is obtained by integrating 

the second-order differential equation for Aiz) of an observing beam from the source at 

□ □ 

redshift 2 ; to the observer at z = 0 [21|, : 

d 


(i+.m.)- 


{1 + z)'^E{z) — \/^aJ^ + + z)^\/ A{z) — 0, (1) 


where E{z) is the reduced Hubble parameter at redshift ^ 

E{z) = ^ = (1 -t- z)^/l + VtmZ + + z)~‘^ — 1], 


( 2 ) 


and the phenomenological parameter rj = 1 — is the so-called dumpiness or smooth¬ 
ness parameter which quantihes the amount of matter in clumps relative to the amount 
of matter uniformly distributed. The required boundary conditions for Eq. ([T]) are 


\/~A |2=o— 0, 


|2=0= -V^- 


dz 

where 6fl is the solid angle of the beam. By using an approximate change of variables 


(3) 


h{A, z) = (1 -F z) 



Ciz) = 


Qr 


1 - 


sn' 


(4) 

(5) 


Eq. ([T]) can be transformed into a hypergeometric equation 

y-^.d^h A 7 .dh +1) 

(i-C)C5^ + (5-gC)5^ + ^^ = o, 


36 


( 6 ) 


The resulting luminosity distance is then given by 


diiz) = (1 + 2;)h(C(0)). 


(7) 
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Expressed in terms of hypergeometric functions, Eq. ([7]) becomes 


dbiz^Qm,^) — —57^-^—[1 + r2m^(3 + 3z + 

Qrn (1 + 2z/) 


X 


2Fi 


V 3 — z/ 5 — 2z/ 


1 - 


6 ’ 6 ’ 6 ' I + VLmz{‘i + z"^) 


X 2F1 


l + v 4 + z/ 7 + 2z/ 


6 


6 


6 


1-a 


- [1 + n^z{?> + 3z + z2)]-d+2-)/6 2 F 1 


z/ 3 — z/ 5 — 2z/ 


6 


6 


1 -a 


X 2^^! 


1 + z/ 4 + Z/ 7 + 2z/ 


1 - 


6 ’ 6 

The parameter v presented in Eqs. (0 
by 


( 8 ) 


’ 6 ' I + VLmZ{?, + ?,Z + z"^) 

and E] corresponds to the dumpiness parameter rj 


ri = x(3 + z/)(2-z/). 
0 


(9) 


The range for z/ is 0 < z/ < 2, where z/ = 0(7 = 1) is related to a FLRW fluid, while 
z/ = 2(7 = 0 ) to a totally clumped case. 

Actually, the ZKDR approach has been criticized by several authors (e.g., a few de¬ 


tailed comments gathered in Ref. 


33|). However, so far, confrontations of the ZKDR 


luminosity distance with observations have not led to conclusive results in the sense of 
totally excluding this model. Moreover, we should keep in mind that most previous tests 
in this held were somewhat dependent on the standard dark energy model (the hat ACDM 
or tcCDM). Therefore, it is necessary to clarify the validity and the scope of the ZKDR 
luminosity distance in describing the Universe in a model-unbiased way. Here, we follow 
the simplest treatment, where rj is assumed to be a constant. 


III. SAMPLES AND RESULTS 


We carry out analysis by using the latest observations of standard candles, including 


the joint light-curve analysis of the SDSS-H and SNLS supernova samples 


48j-which is 


referred to as JLA SN la in the literature-and the long gamma-ray bursts reported in 


Ref. 


49l |. Descriptions for the samples, methodology, and results are presented in this 


section. 
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A. Type la supernovae 


The cosmic acceleration was discovered 16 years ago by measuring accurate distances 


to distant SNe la 


51 


53l |. The reason for the acceleration remains uncertain and a large 


experimental effort in observational cosmology has been driven to reveal the mechanism of 
this ostensibly counterintuitive phenomenon. By precisely mapping the distance-redshift 
relation up to redshift 2 ; ~ 1, SNe la remain, at this stage, the most promising probe 
of the late-time history of the Universe. Because of the variability of the large spectra 
features, distance estimation for SNe la is based on the empirical observation that these 
events form a homogeneous class whose remaining variability is reasonably well captured 
by two parameters 5^. One of them characterizes the stretching of the light curve {Xi 
in what follows), and the other describes the color at maximum brightness (C in what 
follows). 

With the assumption that SNe la at all redshifts with the identical color, shape and 
galactic environment have, on average, the same intrinsic luminosity, the distance esti¬ 
mator (distance modulus: /i = 5 log + 25) used in most cosmological analysis is 
quantihed by a linear model. 


Pb(«, /d; M) = mg — M + a X Xi — (3 x C, 


( 10 ) 


where mg is the observed peak magnitude in the rest-frame B band, and a and /3 are 
nuisance parameters which characterize the stretch-luminosity and color-luminosity re¬ 
lationships, corresponding to the well-known broader-brighter and bluer-brighter rela¬ 
tionships, respectively. The value of M is another nuisance parameter representing the 
absolute magnitude of a hducial SNe la. In general, a and (3 are left as free parameters 
(on the same weight as cosmological parameters) that are determined in the global £t in 
the context of standard dark energy scenario to construct the Hubble diagram for SNe 
la. It should be noted that cosmological implications derived from this Hubble diagram 
for other nonstandard models, which are different from the standard ACDM (or wCDM) 
scenario used to carry out the global £t, are model biased. 

In order to achieve model-unbiased constraints on the dumpiness of the Universe, we 
should £t the light-curve htting parameters {a and (3) and the model parameters (fl^ and 
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z/) simultaneously to construct a Hubble diagram of SNe la in an inhomogeneity-allowed 
scenario by confronting the distances estimated from SNe la observations via Eq. oa 
with the ones predicted from the ZKDR luminosity distance model, 

/Umod(^;6'l,/Uo) = 51ogio[E>L(^;6'l)] + ho- (11) 


Here -Dl is the Hubble-constant free luminosity distance, O-^ represents the model param¬ 
eter vector (Hm, and /xq = 5 logio[c/-Ho] + 25. For the latest JLA SN la, the standard 
function is given by 


2/ n n \ \ '' [hmod(-^ij ^1) ho) hB,j(^ 2 ) Af)] 


x^(ho; Af; 6 ^ 1 , 6 ^ 2 ) — ^2 


CT. 


( 12 ) 

i=i "Md 

where 62 denotes the vector of light-curve htting parameters (a,/5) and is the error 
on the distance modulus for the ith SNe la. It should be noted that we take only the 
statistical uncertainties into account and they are also dependent on the light-curve fitting 
parameters. In order to marginalize over the nuisance parameters, Hq and M, we expand 
the function with respect to /To = ho + Af as |55l-l57| 


X^(6'i, 02 ] ho) — A — 2fioB + 


(13) 


where 


740 


4(»i,»2) = 5^ 


i=l 

740 


s( 9 i,e 2 ) = 5 ^ 


2=1 

740 


[/^mod('^2 7 ^ 1 ; /^O 0 ) /^B,2(^2 7 0 )] 

[/^mod ('^27 ^ 1 : /^O 0 ) 2(^21 ^ 0 )] 


C(»2) = 5; 


2 
=1 


Equation flT3l) has a minimum a.t fio = B/C, and it is 

R 2 

x\0i,02)=A-—. 


( 14 ) 

( 15 ) 

( 16 ) 


(17) 


Therefore, we can minimize x^{0i, 62 ) to get rid of the dependence on nuisance parameters. 

The constraint on the light-curve htting parameters vector is presented in Fig. [H The 
best ht value is (a, /3) = (0.13, 3.17), which is marginally compatible with the result 








estimated in the flat ACDM at a la confldence level. By applying a minimization of 
X^, we can get an estimation for /Tq which is a combination of Hq and M. Here, we 
break the degeneracy by fixing Hq = 70 km s~^ Mpc“^ and obtain M = —19.08. With 
the constraint on {a, /3) and estimation of M, an indicative Hubble diagram in the 
framework of the ZKDR luminosity distance model is constructed and shown in Fig. O 
Moreover, results for confldence regions constrained in the {flm, plane are presented 
in Fig. |3]and Tab. Ull We suggest that the dumpiness parameter rj is poorly constrained, 
being bounded on the interval 0.16 < r] < 1.00 within a la confldence level. However, a 
tighter constraint is obtained for the matter density parameter firm being restricted on the 
interval 0.25 < flm < 0.37(lcr). These are very similar to what was obtained in previous 
analyses [29|, l32|, but quite different from the results included in Ref. [33|. That is, our 
unbiased tests slightly indicate an inhomogeneity and the standard FLRW cosmology is 
consistent with SNe la observations within a Icr confldence level. 



a : (0.13, 3.17) Clumpy universe 


★ : (0.14, 3.14) ACDM 
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FIG. 1: Constraints on the light-cnrve fitting parameters, a and /3, from the global fit in the 
context of a clumpy Universe. The triangle and star represent the best fits when the ZKDR 
approximation and the standard ACDM framework are considered, respectively. 
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FIG. 2: Hubble diagram of the standard candles constructed from the global fit in the context of 
a clumpy Universe. The distance modulus redshift relation of the best-fit ZKDR approximation 
for a fixed Hq = 70 km s“^ Mpc“^ is shown as the solid line. 


B. Long gamma-ray bursts 


Gamma-ray bursts (GRBs), which are the most intensive explosions observed in the 
Universe and thus are visible across much larger distances than SNe la, are deemed as 
a potential probe to explore the Universe at higher redshift, a redshift of at least 6 and 


up to even z = 10 [58 


6l| . Specihcally, relations between the luminosity/energy and 


the measurable properties of the prompt gamma-ray emission imply that GRBs may be 
appropriate candidates for cosmological standard candles. In the past few years, several 
empirical luminosity relations have been statistically inferred from observations. For 
instance, several two-variable relations: the relation between spectral lag an d luminosity 
{'Tiag — L) 62|], the relation between variability and luminosity (V — L) 63|,l6^, the relation 
between peak spectral energy and luminosity (Epeak—L) 65|,|^, the relation between peak 
spectral energy and collimation-corrected energy (F^peak — E^) |^, the relation between 
the minimum raising time in the GRB light curve and luminosity (trt ~ Fl 68l|. and the 
relation between peak spectral energy and isotropic energy (F^peak — F^ 7 ,iso) |69|-have been 
successfully deduced from observations. Meanwhile, a few multivariable relations have 
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FIG. 3: Confidence regions in the (Hmj plane for the model with a ZKDR luminosity distance 
constrained from the JLA SN la. 

also been obtained, such as the connection between Eiso, i^peak, and the break time of the 


optical afterglow light curves (tb) l2Q|, the correlation between the luminosity, Ep^nk, and 
the rest-frame “high-signal” time scale (To, 45 ) TlJ. Moreover, these luminosity relations 


have been proposed to calibrate GRBs as distance indicators (see, e.g.. Refs, 
reviews). 


68 


721 for 


In particular, in Refs. 32|, |33|, distances of GRBs used to constrain the dumpiness of 


the Universe are obtained by calibrating their luminosity relations with low-redshift SNe 
la 


39 


41 


46| . However, it is necessary to make clear that distances of SNe la quoted to 


calibrate luminosity relations are estimated from a global ht in the frame of a standarc 
dark energy model. In other words, the distances of GRBs given in Refs. 


39 


41 


46 


are still somewhat dependent on the standard dark energy model and thus subsequent 
tests for the inhomogeneity of the Universe derived from them are model biased. In this 
work, we construct the Hubble diagram of 116 long GRBs 0 in the framework of an 
inhomogeneous Universe by calibrating their luminosity/energy relations in the global ht 
where the context of the ZKDR luminosity distance model is considered. This Hubble 
diagram can then lead to an unbiased examination of the dumpiness of the Universe. In 
Ref. [49|, six luminosity correlations (riag — L, V — L, Ep^ak — L, Epe^k — E-y, trt — L, 
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Epeak — -E' 7 ,iso) have been derived from the latest observations of 116 long GRBs. In their 
work, it was also found that the intrinsic scatter of the V — L correlation was too large 
to infer an inherent correlation between these two quantities using the currently observed 


GRB events. What is more, the luminosity correlations i^peak — and Ep^ak 


E ■ 

-^7,ISO 


mirror almost the same physics, we should include one of them to avoid strong correlation 
among the luminosity correlations. Therefore, we choose the Epaak — E.y correlation, which 
has a smaller intrinsic scatter, and then use the rest four correlations for the following 
analysis. The same as previous works that derived cosmological implications from GRBs, 
we use only the subsample at z > lA for the complimentary redshift range to the SN la. 

The remaining four luminosity correlations involved in this paper are 


log 


1 erg s 


— = ai + bi log 


’ Rag(l + Z ) 
0.1 s 


-1 


(18) 


log 


1 erg s 


— = 02 + 62 log 


.hipeak(l T 

300 keV 


(19) 


log 


E^ 


1 erg 


= 03 + 63 log 


EpeakiX T 
300 keV 


( 20 ) 


log 


— = O4 + 64 log 


'rRT(l + 2:) ^ 

OT^ 


( 21 ) 


1 erg s 

where o and b are the intercept and the slope of the relation, respectively. In these 
correlations, the isotropic peak luminosity L is given by 


L = AttSR 


bolo? 


( 22 ) 


where Pboio is the bolometric flux of gamma rays in the burst. The isotropic energy 
released in a burst is 

^ 7 ,iso = 47rci|S'boio(l + (23) 

where S'boio is the bolometric fluence of gamma rays in the burst at redshift ; 2 . The total 
collimation-corrected energy can be calculated by 


E^ £'.yjso(l COS^jet)) 


(24) 
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where ^jet is the opening angle of the jet. 

In order to completely avoid any circularity and obtain model-unbiased constraints on 


65 


, we separately calibrate each lumi- 


the dumpiness of the Universe from GRBs 
nosity relation, Eqs. by carrying out a similar simultaneous global htting route 


presented in the above subsection. Results are shown in Tab. [H Here, cxint is the system¬ 
atic error and it can be estimated by Ending the value such that an fif fo each relation 
calibration curve produces a value of reduced of unity j68|. This quantity accounts 


the extra scatter of the luminosity relations. In this global fitting route, we marginalize 
the nuisance parameter Hubble constant by fixing Hq = 70 km s“^ Mpc“^. Following the 
method about uncertainty calculation and distance estimation from calibrated luminosity 


relations 


46 


68|, as shown in Fig. [2l we construct a Hubble diagram of GRBs in the 


context of the ZKDR luminosity distance scenario. In addition, results concerning the 
constraints on model parameters are presented in Fig. [Hand Tab. [TTl It is suggested that a 
Universe composed only by homogeneously distributed matter is strongly favored by GRB 
observations. This is greatly different from what was obtained in previous works j32l. l33|. 

Finally, we perform a joint analysis from the combination of JLA SN la and long GRBs. 
Results are displayed in Fig.[5]and Tab.HTl Within a la confidence level, the matter density 
parameter is restricted in the interval 0.32 < Qm < 0.36 and the smoothness parameter 
is bounded in the interval 0.98 < rj < 1.00. It is shown that the constraint on the matter 
density parameter is mainly dependent on SNe la observations while the estimation of 
the smoothness parameter is basically determined by the long GRBs. The fact that high 
redshift GRBs prefer a homogeneous Universe with a great significance of probability can 
be understood as follows: they explore much larger scales of the Universe and should 
contribute to diminishing the corresponding space parameter. That is, since the Universe 
is more homogeneous on larger scales (a higher redshift), higher value of the smoothness 
parameter rj is favored. In addition, it should be noted that, although large redshift GRBs 
are very important for the tests of the dumpiness parameter, there are only four GRBs 
at redshift larger than 5. 
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Luminosity relation 

a(lcj) 

6(lcj) 

Rnt 

^(zgrb > 1-4) 

Rag -h 

52.60 ± 0.04 

-0.76 ±0.06 

0.12 

26 

Epeak L 

52.10 ±0.04 

1.38 ±0.12 

0.16 

62 

-E^peak -E '7 

50.36 ±0.07 

1.56 ±0.20 

0.01 

12 

trt — L 

52.95 ±0.05 

-1.03 ±0.13 

0.16 

36 


TABLE I: Summary of the constraints on luminosity relations of GRBs from the global fit in 
the context of a clumpy Universe. 
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FIG. 4: Confidence regions in the (14^, plane for the model with a ZKDR luminosity distance 
constrained from the long GRBs. 


Sample 


i/(lo-) 

r]{la) 

JLA SN la 

0.25 <nm< 0.37 

0.00 < 1/ < 1.80 

0.16 <7] < 1.00 

Long GRBs 

0.38 <nm< 0.49 

0.00 < 1/ < 0.48 

0.88 <!]< 1.00 

Joint analysis 

0.32 < Dm < 0.36 

0.00 < 1/ < 0.12 

0.98 < 7? < 1.00 


TABLE 11: Summary of the unbiased constraints on model parameters in the ZKDR luminosity 
distance from observations of standard candles. 
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FIG. 5: Confidence regions in the (Hmj plane for the model with a ZKDR luminosity distance 
constrained from the combination of JLA SN la and long GRBs. 

IV. CONCLUSIONS AND DISCUSSIONS 


In the era of precision cosmology, where one aims at determining the cosmological pa¬ 
rameters at the percent level, distance estimations for standard candles and rulers with 
increasing accuracy are expected to provide powerful constraints on dark energy or other 
fundamental dynamical parameters. However, it is necessary to be aware of the physical 
hypothesis underlying these probes when we proceed with such a program. As far as 
we know, the Universe is effectively inhomogeneous at least in the small-scale domain. 
Furthermore, notice that even the large-scale homogeneity also has been challenged 0 . 
In this topic, the method based on the ZKDR luminosity distance is a simple alternative 
and is usually applied to quantitatively assessing the influences of the dumpiness on the 
light propagation. In the past few years, there has been a rich literature concer ning the 


constraints on the smoothness parameter from observations of standard candles 30 


32- 


34 l5nl |. However, we should keep in mind that distances of SNe la applied to test the 


inhomogeneity were estimated from the global £t in the context of a standard homo¬ 
geneous dark energy model, i.e., the flat ACDM or tcCDM model. Therefore, in these 
previous analyses, constraints on the smoothness parameter from the distance modulus 
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of SNe la were somewhat model biased. Meanwhile, results obtained from GRBs suffered 
the same problem since the distances of them were determined by calibrating luminosity 
relations with low-redshift SNe la. 

In this paper, we hrst construct Hubble diagrams for SNe la and GRBs by calibrating 
the light-curve htting parameters and luminosity relations, respectively, in the global 
£t where the context of the ZKDR luminosity distance model is taken into account. 
And then, these Hubble diagrams can lead to unbiased tests for the inhomogeneity of 
the Universe. For the JLA SN la, as shown in Fig. [3l constraint on the smoothness 
parameter is not stringent and slightly implies a locally inhomogeneous background, while 
the matter density parameter is well constrained, being bounded in the interval 0.25 < 
VLm < 0.37(1(t). For the long GRBs, as shown in Fig. 01 the Universe with matter 
uniformly distributed is favored with a high conhdence level. This is completely different 
from what was obtained in Ref. 33|]. Finally, we perform a joint analysis which provides 


good constraints on both model parameters. At a la conhdence level, the intervals are 
0.32 < flm < 0.36 and 0.98 < rj < 1.00. It is suggested that the constraint on the 
matter density parameter is mainly based on the observations of low-redshift SNe la, 
while the test for the dumpiness parameter is primarily determined from the observations 
of high-redshift GRBs. Just as expected, the investigation on the inhomogeneity was 
very sensitive to the scales explored by the observations, i.e., the Universe should be 
more homogeneous on larger scales. These also may be an indication that the ZKDR 
approximation remains to be a precise description for the luminosity distance-redshift 
relation in a locally inhomogeneous Universe with the cosmological constant. 

Frankly, it should be pointed out that constraints on the model parameters from low- 
redshift SNe la and high-redshift GRBs are somewhat inconsistent. This inconsistency 
may imply that the assumption with the smoothness parameter rj being a constant is 
not accurate enough to £t the practical observations. That is, the smoothness parameter 
7 ] might evolve with cosmic time (or redshift). Moreover, the intrinsic scatters in GRB 
observations may also lead to this tension. Therefore, in the near future, a more precise 
and larger sample of high-redshift GRB data (even some other distance measurements 
with new methods, e.g., extremely luminous active galactic nuclei readily observed over 
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a range of distances from ~ 10 Mpc to z > 7 j74l476l|i and a plansible extension of the 


ZKDR approach are expected to perform more accnrate tests for the inhomogeneity and 
contribntion of matter in the Universe. 
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